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Computational  fluid  dynamics  calculations  using  high-performance  parallel  computing 
were  conducted  to  simulate  the  prestall  flow  of  a  transonic  compressor  stage,  NASA 
compressor  Stage  35.  The  simulations  were  run  with  a  full-annulus  grid  that  models  the 
3D,  viscous,  unsteady  blade  row  interaction  without  the  need  for  an  artificial  inlet  dis¬ 
tortion  to  induce  stall.  The  simulation  demonstrates  the  development  of  the  rotating  stall 
from  the  growth  of  instabilities.  Pressure  rise  performance  and  pressure  traces  are  com¬ 
pared  with  published  experimental  data  before  the  study  of  flow  evolution  prior  to  the 
rotating  stall.  Spatial  fast  Fourier  transform  analysis  of  the  flow  indicates  a  rotating 
long-length  disturbance  of  one  rotor  circumference,  which  is  followed  by  a  spike-type 
breakdown.  The  analysis  also  links  the  long-length  wave  disturbance  with  the  initiation  of 
the  spike  inception.  The  spike  instabilities  occur  when  the  trajectory  of  the  tip  clearance 
flow  becomes  perpendicular  to  the  axial  direction.  When  approaching  stall,  the  passage 
shock  changes  from  a  single  oblique  shock  to  a  dual  shock,  which  distorts  the  perpen¬ 
dicular  trajectory  of  the  tip  clearance  vortex  but  shows  no  evidence  of  flow  separation 
that  may  contribute  to  stall.  [DOI:  10.1115/1.2812968] 


Introduction 

Advances  in  gas  turbine  engine  technology  have  made  possible 
the  production  of  highly  efficient  engines.  However,  the  useful 
operating  range  of  a  gas  turbine  engine  is  limited  by  the  stable 
operation  of  the  compressor.  The  limiting  instability  can  be  char¬ 
acterized  by  either  rotating  stall  or  surge.  Surge  is  the  zeroth  order 
planar  oscillation  mode  resulting  from  overall  imbalance  of  the 
entire  propulsion  system,  thus,  analysis  would  have  to  include  full 
volume  and  ducting.  This  violent  system  instability  may  be  pre¬ 
ceded  by  a  rotating  stall,  which  sets  up  alternating  stress  on  the 
blades.  Rotating  stall  and  surge  can  cause  destructive  damage  to 
engines  and  should  be  avoided.  Understanding  the  physical 
process  of  the  onset  of  rotating  stall  can  aid  in  the  design  of  a 
better  control  method  to  stabilize  the  system  and  enhance  engine 
performance. 

Stall  inception  has  been  an  active  research  effort  for  the  past 
five  decades.  The  events  leading  to  rotating  stall  have  been  tradi¬ 
tionally  classified  according  to  two  different  types  of  wave  distur¬ 
bances  moving  around  the  annulus:  long-length  (modal)  and 
short-length  (spike)  waves.  Long-length  wave  disturbances  are 
those  with  length  on  the  order  of  the  rotor  circumference.  This 
type  of  disturbance  is  caused  by  the  characteristic  frequencies  of 
the  compression  system.  Theories  to  model  the  long-length  insta¬ 
bility  were  developed  by  Moore  and  Greitzer  [1].  The  assumption 
is  that  the  disturbance  length  scale  is  much  larger  than  one  blade 
pitch,  and  that  the  entire  blade  row  can  be  modeled,  both  spatially 
and  temporally,  without  the  knowledge  of  the  detailed  flow  inside 
the  blade  passages.  This  approach  predicts  the  existence  of  long- 
length  waves  leading  to  stall.  McDougall  [2]  demonstrated  the 
existence  of  the  modal  waves  through  experiments,  and  Gamier 
et  al.  [3]  supported  the  validity  of  the  Moore  and  Greitzer  model 
with  detailed  measurements.  Unsteady  pressure  or  axial  velocity 
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measurements  from  probes  placed  around  the  annulus  are  usually 
used  to  track  the  development  of  modal  waves.  Modal  oscillations 
often  have  an  axial  extent  from  inlet  to  exit.  It  is  typically  a  2D 
( 0,x )  phenomenon,  which  is  not  an  early  form  of  stall  cell  but 
instead  represents  harmonic  oscillations  of  the  flow  field.  It  usu¬ 
ally  appears  when  the  shape  of  the  total-to- static  pressure  rise 
characteristic  becomes  horizontal. 

The  second  type  of  stall  inception  mechanism  is  the  short- 
length,  or  spike,  disturbances.  These  are  disturbances  with  length 
scale  on  the  order  of  blade  passage  breadth.  Spikes  can  be  viewed 
as  embryonic  stall  cells  with  flow  breakdown  in  local  regions 
(Camp  and  Day  [4]);  it  is  a  3D  phenomenon  whose  development 
depends  on  the  flow  structure  within  the  blade  passage.  A  com¬ 
mon  scenario  for  the  spike  stall  inception  process  is  that  localized, 
3D  pockets  of  instability  in  multiple  blade  passages  empt  first  to 
form  small-sized  multicell  rotating  stall.  These  stall  cells  rotates  at 
30-85%  rotor  speed  initially  and  then  merge  into  a  single  rotating 
stall  cell  within  a  few  rotor  revolutions  (Day  [5]).  The  final  speed 
of  the  consolidated  stall  cell  is  nearly  half  or  less  of  the  rotor 
speed  (Cumpsty  [6]).  A  spike  disturbance  more  frequently  occurs 
when  the  slope  of  the  pressure  rise  characteristic  is  still  negative. 
That  is,  the  pressure  still  rises  as  the  mass  flow  reduces. 

The  above  discussions  indicate  that  the  slope  of  the  pressure 
rise  characteristic  may  correlate  with  the  stall  inception  type. 
Spike  and  mode  can  also  be  the  different  aspects  of  the  stall  in¬ 
ception  process  and  can  exist  in  the  same  machine  (Day  [5]).  Day 
indicated  that  modal  waves  may  be  measurable  in  the  machine  for 
as  many  as  200  rotor  revolutions  before  stall.  A  spike  disturbance, 
on  the  other  hand,  once  formed,  can  lead  to  rotating  stall  quickly, 
usually  within  a  few  rotor  revolutions. 

Numerical  simulation  using  computational  fluid  dynamics 
(CFD)  techniques  has  been  shown  as  an  alternative  way  to  ana¬ 
lyze  aerodynamic  problems  and  is  a  complement  to  experiments 
in  the  understanding  of  complex  flow  phenomena.  Flow  structure 
can  be  examined  in  detail  where  experimental  measurements  are 
difficult  or  impossible.  During  the  past  ten  years,  steady-flow 
simulations  have  provided  an  increasingly  accurate  prediction  of 
the  flow  up  to  the  point  of  compressor  stall.  Attempts  to  study  stall 
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through  unsteady  simulations  of  a  periodic  subset  of  the  blades  in 
a  compressor  blade  row  (Hoying  et  al.  [7],  Vo  et  al.  [8],  Davis  and 
Yao  [9])  or  through  reduced-order  unsteady-flow  models  (Gong 
et  al.  [10])  have  reached  valuable  findings  of  the  flow  structure 
leading  to  stall.  However,  since  the  temporal  flow  field  variations 
that  occur  during  stall  inception  are  not  harmonics  of  blade  pass¬ 
ing  frequency,  a  full-annulus  model  should  be  used  to  remove  the 
constraint  of  the  periodic  model  in  order  to  study  the  transition 
from  a  steady  flow  into  the  unsteady  stalling  flow,  and  to  resolve 
flow  features  of  length  scale  on  the  order  of  the  rotor  circumfer¬ 
ence.  Such  simulations  have  been  done  two  dimensionally  (He 
[11],  Saxer-Felici  et  al.  [12]),  but  stall  and  the  instabilities  leading 
to  stall  are  inherently  three-dimensional  phenomena.  Predicting 
stall  phenomena  may  thus  require  three-dimensional  unsteady 
full-annulus  simulations,  which  is  a  daunting  computational  en¬ 
deavor  requiring  considerable  computational  resources.  Due  to  the 
rapid  advance  of  high-performance  computing  in  the  past  decade, 
advanced  CFD  codes  have  benefited  from  the  almost-unlimited 
computing  potential  of  parallel  computing.  Recently,  Chen  et  al. 

[13]  demonstrated  the  ability  to  conduct  a  full-annulus  3D  un¬ 
steady  computation  of  an  entire  compressor  stage  while  Hah  et  al. 

[14]  also  conducted  a  full-annulus  3D  computation  of  an  isolated 
rotor.  This  trend  indicates  that  large-scale  simulation  using  paral¬ 
lel  computing  is  within  reach. 

Previous  comparisons  of  the  time  averages  of  unsteady  phase- 
lag  and  partial  annulus  predictions  to  measurements  have  shown 
the  validity  of  the  turbo  code  to  capture  the  ID  and  2D  charac¬ 
teristics  of  this  compressor  (Hathaway  et  al.  [15]).  The  intent  of 
this  work  is  to  use  the  available  experimental  data  to  provide 
assessment  of  the  code  and  to  show  evidence  that  the  computa¬ 
tional  model  reasonably  captures  the  inherent  behavior  of  com¬ 
pressors  as  they  are  throttled  into  stall.  It  is  recognized  that  some 
features  of  the  actual  tested  compressor  are  not  captured  by  the 
computational  model,  for  example,  the  leakage  gap  at  the  rotor/ 
stator  hub  interface,  the  actual  blade  tip  profile,  and  the  geometry 
of  the  upstream  and  downstream  ducting  between  the  plenum  and 
throttle  valve  locations.  Although  these  features  might  influence 
the  actual  compressor  stall,  they  are  not  considered  requisite  to 
understanding  fundamental  stall  inception  characteristics. 

Axial  Compressor,  Stage  35 

NASA  research  single-stage  compressor  Stage  35,  representa¬ 
tive  of  the  advanced  transonic  core  compressor,  is  used  in  this 
investigation.  A  schematic  of  the  compressor  test  rig  and  compres¬ 
sor  geometry  is  provided  in  Fig.  1  and  the  design  parameters  are 
provided  in  Table  1.  Stage  35  produces  1.8  total  pressure  ratio  at 
a  mass  flow  rate  of  20.2  kg/ s  at  the  design  speed  of  17,189  rpm. 
Details  of  the  Stage  35  geometry,  operating  conditions,  and  the 
experimental  results  are  provided  by  Reid  and  Moore  [16]. 

Relevant  Experimental  Work 

Time-resolved  pressure  measurement  of  Stage  35  near  stall  was 
reported  by  Bright  et  al.  [17].  The  experiment  was  performed  at 
85%  design  speed.  Several  configurations,  including  a  clean-inlet 
case,  were  reported  in  that  paper.  For  the  interest  of  this  study  we 
will  focus  on  the  clean-inlet  case.  The  compressor  was  transi¬ 
tioned  into  stall  using  a  continuous  throttle  maneuver.  Time- 
resolved  pressure  was  recorded  by  high-frequency-response  trans¬ 
ducers  at  eight  locations  equally  spaced  around  the  annulus. 
Figure  2  reproduces  the  pressure  traces  of  the  clean-inlet  case  for 
180  revolutions  prior  to  stall  reported  in  that  paper  (the  original 
data  are  not  available  to  us).  Based  on  the  pressure  disturbance, 
Bright  et  al.  [17]  commented  that  the  compressor  with  the  clean- 
inlet  stalls  due  to  modal  instabilities,  and  at  40  revolutions  before 
stall,  there  appear  to  be  spikes  on  top  of  the  modal  wave.  They 
concluded  that,  “This  indicates  that  our  high-speed  stage  under 
normal  stalling  condition  is  considered  a  ‘modal’  machine,  since 
modal  instabilities  dominate  the  pre- stall  behavior  with  some 
added  pip  instabilities  present  just  before  stall.” 
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Fig.  1  NASA  high-speed  compressor  Stage  35 


In  the  following,  we  take  a  closer  look  at  these  experimental 
data.  The  shapes  of  the  modal  disturbances  indicated  in  Fig.  2 
resemble  sinusoidal  waves  with  a  temporal  period  of  ten  rotor 
revolutions.  The  speed  of  the  modal  disturbances,  as  indicated  by 
the  slope  of  Line  A,  connecting  the  peaks  of  the  pressure  wave,  is 
zero.  This  indicates  that  the  modal  disturbance  is  caused  by  axi- 
symmetric  longitudinal  system  disturbances.  The  system  includes 
all  ducting  and  volumes  of  the  experimental  facility.  This  modal 
disturbance  is  the  surge  mode  that  oscillates  at  a  frequency  of  10% 
rotor  speed — equivalent  to  the  ten  rotor  revolutions  period  it  cov¬ 
ers.  Unlike  a  typical  rotating  modal  wave,  the  surge  mode  does 
not  rotate.  The  surge  mode  is  a  compression  system  wide  reso¬ 
nance  that  cannot  be  modeled  with  the  current  simulation  because 
the  computational  grid  does  not  include  the  ducting  of  the  test 
facility,  which  is  quite  substantial,  as  shown  in  Fig.  1. 

As  the  compressor  throttles  closer  to  stall  at  about  ten  revolu¬ 
tions  before  stall,  there  is  a  secondary  disturbance,  much  like  a 
spike  in  the  time  scale  of  Fig.  2.  Looking  more  closely,  this  sec¬ 
ondary  disturbance  covers  a  circumferential  distance  of  one  rotor 
circumference,  as  indicated  in  the  enlarged  image.  It  sits  on  top  of 


Table  1  Design  parameters  of  Stage  35 


Rotor  rpm  at  100%  speed 

17,188.7 

Tip  speed  (m/s) 

454.456 

Hub/tip  radius 

0.7 

Rotor  aspect  ratio 

1.19 

Stator  aspect  ratio 

1.26 

Number  of  rotor  blades 

36 

Number  of  stator  blades 

46 
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Pressure  Traces  for  No  Control,  No  Blowing  Stage  35 


1 1 - • - * - i - i - \ - « - 

-140  -100  -60  -20  0 

Rotor  Revolutions 

Fig.  2  Experiment  pressure  traces  for  clean  inlet,  85%  design 
speed,  Stage  35,  Bright  et  al.  [17] 


the  surge  mode  and  rotates  around  the  annulus  at  about  80%  rotor 
speed,  indicated  by  Line  B.  Because  the  length  scale  of  this  dis¬ 
turbance  is  the  entire  annulus,  which  is  much  larger  than  that  of  a 
typical  spike,  it  is  classified  as  a  mode.  For  discussions  in  the  rest 
of  the  paper,  this  rotating  long-length  disturbance  is  also  referred 
to  as  a  rotating  modal  disturbance.  This  long-length  disturbance  is 
followed  by  the  eventual  flow  breakdown.  Due  to  the  picture  qual¬ 
ity,  we  are  unable  to  discern  the  transition  between  the  start  of  the 
rotating  long-length  disturbance  and  the  final  breakdown.  The 
long-length  disturbance  has  a  length  scale  of  one  rotor  circumfer¬ 
ence.  The  rotating  long-length  disturbance,  because  of  its  length 
scale,  can  only  be  resolved  by  a  full-annulus  simulation. 

The  simulation  reported  in  this  work  is  part  of  an  overall  effort 
to  examine  steady  tip  injection  stall  flow  control  for  Stage  35  at 
design  speed  with  the  clean  inlet.  As  such,  the  simulation  was 
conducted  at  design  speed.  To  our  knowledge,  there  are  no  pres¬ 
sure  trace  data  available  at  design  speed  as  that  shown  in  Fig.  2. 
Therefore,  we  will  use  the  85%  speed  case  to  compare  with  the 
simulated  results.  While  the  speeds  are  different,  the  two  cases  are 
both  transonic  compressors  and  share  common  features  as  we  ex¬ 
plore  the  flow  behavior  prior  to  stall.  The  experimental  character¬ 
istics  at  85%  and  100%  speeds  are  compared  in  Fig.  3.  Both  cases 
stall  at  the  peak  of  static  pressure  rise.  More  importantly  is  that 
their  slopes  near  stall  are  similar.  These  similarities  justify  the  use 
of  the  85%  speed  case  for  comparison. 

Numerical  Methodology 

turbo  is  a  physics-based  simulation  tool  for  multistage  turbo¬ 
machinery.  The  solver  computes  the  fluid  conservation  laws  with¬ 
out  ad  hoc  modeling  of  any  flow  phenomenon  other  than  models 
required  for  turbulence.  This  code  solves  the  unsteady  Reynolds- 
averaged  Navier-Stokes  equations  and  a  decoupled  k-s  turbulence 
model.  The  code  is  implemented  in  a  portable,  scalable  form  for 
distributed-memory  parallel  computers  using  message  passing  in¬ 
terface  (MPI).  The  parallel  implementation  employs  domain  de¬ 
composition  and  supports  general  multiblock  grids  with  arbitrary 
grid-block  connectivity.  The  solution  algorithm  is  a  Newton  itera¬ 
tive  implicit  time-accurate  scheme  with  characteristics-based 
finite-volume  spatial  discretization.  The  Newton  subiterations  are 
solved  using  a  concurrent  block-Jacobi  symmetric  Gauss-Seidel 
(BJ-SGS)  relaxation  scheme.  Because  all  of  the  fundamental  fluid 
mechanics  are  computed,  the  code  is  capable  of  capturing  the 
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Fig.  3  Comparison  of  characteristics  at  85%,  Bright  et  al.  [17], 
and  100%,  Weigl  et  al.  [27],  speeds,  Stage  35 


nonlinear  characteristics  of  the  flow  fields  of  interest.  With  the 
actual  modeling  of  blade  rows  in  relative  motion,  this  code  is 
capable  of  computing  the  unsteady  interactions  between  blade 
rows.  Details  of  the  flow  solver  can  be  found  from  Chen  and 
Whitfield  [18].  The  approach  to  parallelization  for  large-scale, 
complex  problems  is  discussed  by  Chen  and  Briley  [19].  In  the 
following,  we  discuss  several  issues  applied  to  the  code  in  the 
simulation  of  rotating  stall. 

Full  Compressor  Stage  Simulation.  The  intent  to  simulate  the 
inception  and  maturation  of  long-wavelength  circumferential 
modes,  with  wave  numbers  as  low  as  unity,  requires  full-annulus 
simulations.  A  time-accurate  sliding  interface  (Chen  and  Barter 
[20])  is  implemented,  in  which  pressure  waves  in  all  directions 
(axial,  radial,  and  circumferential)  are  passed  unaltered  between 
the  rotor  and  the  stator.  This  attribute  is  necessary  to  examine  the 
temporal  and  spatial  responses  associated  with  blade  row  interac¬ 
tion  during  stall  inception. 

Inlet  Boundary  Condition.  An  isentropic  inlet  condition 
(Chen  and  Whitfield  [18])  is  applied,  where  the  span  wise  total 
condition  is  prescribed.  This  boundary  condition  preserves  the  up¬ 
stream  total  conditions  at  test  facility  level  but  can  potentially 
reflect  outgoing  pressure  waves  back  to  the  computational 
domain. 

Exit  Boundary  Condition.  The  traditional  choice  of  exit 
boundary  for  stable  flow  is  the  radial  equilibrium  condition.  This 
boundary  condition  uses  a  preset  exit  static  pressure.  In  near-stall 
or  during-stall  conditions,  compressor  pressure  rise  drops  and  a 
preset  exit  pressure  cannot  match  the  pressure  drop.  In  this  study, 
we  relieve  the  situation  by  the  use  of  a  “choked”  throttle  model 
that  specifies  corrected  mass  flow  at  the  exit.  This  boundary  con¬ 
dition  allows  variation  of  exit  static  pressure  to  match  the  com¬ 
pressor  exit  mass  corrected  to  the  exit  total  condition. 

Tip  Clearance  Model.  The  tip  gap  was  modeled  by  the  tip 
clearance  model  of  Kirtley  [21].  This  model  ensures  the  conser¬ 
vation  of  mass  and  momentum  through  the  tip  gap  without  the 
need  of  a  tip  clearance  grid.  A  gridding  strategy  suggested  by  Van 
Zante  et  al.  [22]  is  adopted  in  this  study  where  ten  points  are 
placed  in  the  tip  gap  region.  Although  the  detailed  physics  of  the 
development  of  the  leakage  vortex  are  not  captured,  the  tip  clear¬ 
ance  model  was  considered  adequate  for  predicting  the  first  order 
effects  of  the  leakage  vortex  strength,  extent,  and  direction 
(Chima  [23]). 

Gridding.  This  compressor  stage  consists  of  36  rotor  blades 
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1.6 


Fig.  4  Three-blade-row  grid  model  for  the  Stage  35  simulation 


and  46  stator  blades.  For  the  purpose  of  conducting  simulations 
with  tip  injectors  in  another  effort,  we  built  a  three-blade-row  grid 
as  shown  in  Fig.  4.  The  injector  row  has  12  injectors  (not  acti¬ 
vated  in  this  study)  that  are  equally  spaced  in  the  circumferential 
direction.  The  grid  size  of  the  rotor  is  151(x),  71(r),  and  56(0) 
with  81(v)  and  61(r)  points  on  the  blade.  The  grid  size  of  the 
stator  is  141(v),  71(r),  and  79(0)  with  81(v)  and  53  (r)  points  on 
the  blade.  The  rotor  tip  clearance  is  1%  of  rotor  tip  chord.  The 
grids  are  replicated  around  the  annulus  to  generate  a  full-annulus 
grid  with  a  total  of  67  X  106  grid  points.  The  grid  is  then  parti¬ 
tioned  into  328  blocks.  The  simulations  were  run  on  an  IBM  P4 
cluster  at  NAVO  MSRC  using  328  CPUs  with  24  h  (wall  clock) 
computation  time  for  one  rotor  revolution.  In  the  simulation  of  the 
development  of  rotating  stall,  ten  revolutions  were  needed.  Ide¬ 
ally,  each  operating  point  should  be  run  at  least  ten  revolutions  to 
ensure  convergence,  especially  for  points  close  to  stall.  The  time 
for  each  speedline  would  then  be  0.5  X  106  CPU  hours  (six  oper¬ 
ating  points). 

Constrained  by  the  computational  resources  for  the  full-annulus 
simulations,  very  short  upstream  and  downstream  computational 
domains  were  used  here.  This  choice  prompts  two  caveats:  (1) 
The  exclusion  of  the  test  facility  volume  effect  prevents  the  mod¬ 
eling  of  the  surge  mode.  Any  modal  activity  showing  up  in  the 
simulation  should  be  a  result  of  the  compound  effect  of  the  blades 
(both  rotor  and  stator  if  the  interaction  of  the  two  matters)  and  the 
limited  volume  effect  present  in  the  ducting  of  the  current  grid.  (2) 
Limiting  the  axial  extent  of  the  computational  domain  can  poten¬ 
tially  disrupt  long-length  modes  due  to  the  reflective  inlet  bound¬ 
ary.  Because  pressure  waves  (with  length  X)  decay  e~L,x  through  a 
distance  L,  it  is  preferred  to  place  the  inlet  boundary  away  from 
the  rotor  so  that  pressure  waves,  especially  with  long  wave¬ 
lengths,  decay  sufficiently  before  reaching  the  inlet.  Typically,  the 
modal  instability  shows  up  hundreds  of  revolutions  prior  to  stall. 
Complete  simulation  of  the  evolution  of  such  modal  instability  is 
beyond  the  computational  resources  available  in  this  study.  Our 
intent,  therefore,  is  to  examine  the  development  of  the  modal  dis¬ 
turbances  as  they  first  emerge  before  reaching  the  inlet  rather  than 
the  complete  resolution  of  the  modal  stall  process.  In  this  aspect, 
a  shorter  upstream  domain  (2.5  rotor  chords)  is  justified. 
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Results  and  Discussions 

The  difference  of  the  throttling  processes  in  the  simulation  and 
the  experiment  is  that  in  the  simulation,  the  exit  boundary  condi¬ 
tion  is  changed  in  discrete  steps,  while  in  the  experiment,  the 
compressor  transitioned  to  stall  through  a  continuous  throttling 
maneuver.  In  this  simulation,  a  new  throttle  setting  restarts  from 
the  solution  of  a  previous  throttle  setting.  Due  to  computational 
resources,  each  throttle  setting  was  only  run  for  a  few  revolutions 
before  moving  to  the  next.  Generally,  it  was  sufficient  to  reach 
convergence  within  four  revolutions  while  away  from  stall.  As 
stall  was  approached,  more  (six  or  more)  revolutions  were  re¬ 
quired  to  assure  convergence. 

Stage  Characteristics.  The  computed  and  experimental  char¬ 
acteristics  (speedlines)  are  shown  in  Fig.  5.  The  data  used  in  this 
figure  are  the  ratio  of  exit  static  pressure  to  inlet  total  pressure 
versus  mass  rate  corrected  to  inlet  total  condition.  The  computed 
speedline  is  generated  by  incrementally  closing  the  throttle  setting 
(decreasing  exit  corrected  mass  flow)  from  A  to  F.  The  overall 
shape  of  the  computed  characteristic  before  stall  is  similar  to  the 
experimental  data,  so  is  the  slope  of  the  characteristic  near  the 
peak  pressure  (throttle  Point  D).  Since  the  slope  of  the  compressor 
characteristic  correlates  with  the  type  of  stall  inception,  we  con¬ 
sider  the  simulation  reasonable  to  capture  the  stall  inception 
mechanisms. 

Point  D  in  Fig.  5  is  the  last  stable  throttle  setting.  Points  E  and 
F  both  restart  from  Point  D  and  develop  into  a  rotating  stall.  In 
both  cases,  the  exit  pressure  changes  from  the  peak  at  D  into  a 
transitional  peak  and  then  drops  to  a  much  lower  level  during 
rotating  stall,  as  shown  in  Fig.  6(a).  This  demonstrates  the  throt¬ 
tling  effect  of  the  corrected  mass  exit  boundary  condition  by  al¬ 
lowing  the  exit  pressure  to  vary  as  the  flow  condition  changes. 
The  pressure  ratios  shown  in  Fig.  5  for  E  and  F  are  the  transitional 
peak  pressures  before  entering  rotating  stall.  Point  D  is  the  ap¬ 
proximate  stall  point  from  the  experiment.  The  actual  simulated 
stall  point  should  occur  between  D  and  E.  The  actual  stall  point 
depends  on  many  issues,  such  as  tip  clearance  height,  deviation  of 
real  blade  geometry  from  designed  geometry,  as  well  as  the  type 
of  turbulence  modeling  in  the  numerical  simulation.  These  details 
may  not  be  correctly  represented  in  this  simulation.  Even  though 
the  exact  stall  point  is  important,  for  the  intent  of  this  study,  the 
simulation  is  considered  adequate  to  examine  the  general  physical 
trends  leading  to  stall. 

In  the  following  sections,  we  first  present  the  entire  evolution  of 
the  rotating  stall  as  the  throttle  setting  changes  from  D  to  F.  This 
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Fig.  6  Pressure  and  mass  variation  during  rotating  stall  at 
throttle  Setting  F:  (a)  Pressure  ( b )  mass  flow 


is  followed  by  a  closer  examination  of  the  prestall  flow  develop¬ 
ment  from  D  to  E,  with  E  being  closer  to  the  stable  point  D. 

Evolution  of  the  Rotating  Stall.  This  result  was  reported  by 
Chen  et  al.  [13]  and  is  outlined  here  for  completeness.  The  rotat¬ 
ing  stall  was  captured  by  the  time  variation  of  static  pressure  at 
eight  pressure  transducers  circumferentially  spaced  upstream  of 
the  rotor.  The  probes  are  located  44%  chord  upstream  of  the  rotor 
leading  edge  at  98%  span.  The  circumferential  locations  are 
10  deg,  70  deg,  100  deg,  160  deg,  190  deg,  250  deg,  280  deg, 
and  340  deg.  The  history  of  pressure  variations  are  shown  in  Fig. 
7,  in  which  the  pressure  levels  are  offset  by  the  probe  circumfer¬ 
ential  locations.  Because  the  volume  effect  is  not  modeled  in  the 
simulation,  the  pressure  disturbance  of  Fig.  7  should  be  viewed  as 
a  filtered  version  of  Fig.  2  by  removing  the  surge  mode 
oscillations. 

Small,  periodic  pressure  variations  caused  by  the  rotor  blade 
passing,  as  indicated  by  the  thickness  of  the  curves,  are  present  in 
all  curves  of  Fig.  7.  A  smooth  disturbance  is  visible  at  t=0.1T 
(T=time  for  one  rotor  revolution)  traveling  at  full  rotor  speed. 
This  disturbance,  analyzed  in  further  detail  later,  has  a  length  scale 
of  one  rotor  circumference,  thus  is  classified  as  a  modal  distur¬ 
bance.  The  rotating  modal  disturbance  continues  to  grow  in  mag¬ 
nitude  but  the  speed  slows  to  84%  rotor  speed  in  approximately 
one  revolution.  After  2 T,  it  transforms  into  several  spike  distur¬ 
bances  around  the  annulus,  each  with  a  smaller  length  scale  of  a 


Fig.  7  Time  history  of  static  pressure  variation  at  eight  loca¬ 
tions  around  the  annulus  located  44%  chord  ahead  of  the  rotor 
of  Stage  35,  throttling  Setting  F 


few  blade  passages.  The  spike  disturbances  then  grow  in  magni¬ 
tude  and  slow  down  to  43%  in  four  revolutions.  Inlet  and  exit 
mass  flow  histories  during  the  same  period  are  shown  in  Fig.  6(b) , 
in  which  the  mass  flow  remains  relatively  flat  in  the  first  two 
revolutions  at  18  kg/s  before  it  drops  sharply  and  sustains  large 
variations  of  14%  of  the  peak  mass  flow. 

In  Figs.  8  and  9,  the  entropy  on  a  surface  of  revolution  near  the 
casing  is  shown  at  two  instants  in  time.  The  full-annulus  surface  is 
illustrated  by  two  half-annulus  surfaces.  At  /=3.8T,  Fig.  8,  the 
early  phase  of  the  stall  inception,  three  separate  high  entropy  (red 
color)  regions  can  be  seen,  with  the  axial  extent  covering  from 
one  chord  upstream  to  aft  of  the  stator  leading  edges.  The  radial 
extent  of  the  stall  cell  is  shown  in  the  axial  cut  plane  at  16%  chord 
before  the  rotor,  where  the  high  entropy  region  shows  the  stall  cell 
extending  about  30%  span  from  the  casing.  At  t-  10T,  Fig.  9,  the 
stall  cell  evolves  into  the  final  stage  as  its  rotating  speed  becomes 
stabilized.  The  three-cell  stall  has  merged  into  a  single-cell  stall 
with  the  axial  extent  covering  all  the  way  from  the  inlet  of  the 
computation  domain  to  aft  of  the  stator  trailing  edge.  The  radial 
extent  of  the  stall  cell  has  now  extended  about  70%  span  from 
casing. 

The  evolution  of  the  simulated  rotating  stall  cell  is  similar  to 
the  spike  stall  inception  mechanism  commonly  observed  in  that 
localized,  three-dimensional  pockets  of  instability  in  multiple  ro¬ 
tor  passages  break  up  to  form  the  initial  phase  of  a  rotating  stall 
cell,  usually  with  several  small- sized  stall  cells.  These  small  stall 
cells  move  fast  at  the  beginning  and  then  merge  into  a  single 
rotating  stall  cell  within  a  few  revolutions  before  slowing  down  to 
about  half  of  the  rotor  speed  as  the  rotating  stall  becomes  fully 
developed. 

In  the  previous  section,  we  gave  a  brief  description  of  the  over¬ 
all  evolution  of  the  rotating  stall  using  throttle  Setting  F.  In  the 
following  sections,  we  will  examine  the  process  that  leads  to  the 
initial  flow  breakdown,  i.e.,  the  stall  inception  process.  Restarting 
from  the  stable  solution  of  D,  the  simulation  was  repeated  with 
throttle  Setting  E,  the  mass  flow  history  during  the  transition  from 
D  to  E  is  shown  in  Fig.  10.  The  simulation  shows  that  after  four 
revolutions  from  restart,  the  mass  flow  starts  to  decrease,  leading 
to  the  rotating  stall.  Three  time  stamps  marked  El  (1.77), 
E2  (3.27),  and  E3  (4.37)  will  be  referenced  in  the  following 
discussions. 

Development  of  Long-Length  Disturbance.  Both  experimen¬ 
tal  (Spakovszky  et  al.  [24]  and  Suder  et  al.  [25])  and  computa¬ 
tional  studies  (Chen  et  al.  [13],  Davis  and  Yao  [9])  indicate  that 
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Fig.  8  Entropy  of  Stage  35  during  stall  inception,  f=3.8T:  Left 
half  annulus  and  right  half  annulus,  16%  chord  upstream  of 
rotor  leading  edge 


the  rotating  stall  of  this  compressor  initiates  in  the  rotor  tip  region, 
so  the  analysis  will  focus  on  the  rotor  tip.  Instantaneous  pressure 
profiles  around  the  annulus  near  the  rotor  tip  leading  edge  are 
shown  in  Fig.  11.  The  three  curves  correspond  to  different  phases 
at  throttle  Settings  D  and  E  during  the  transition  from  stable  op¬ 
eration  to  stall  (Fig.  10).  Phase  1  is  the  last  stable  point  before 
stall,  which  will  be  used  as  the  base  line  for  comparison.  Phases  2 
and  3  are  the  same  throttle  Setting  E,  but  at  two  different  instants: 
one  before  stall  (Time  El)  and  the  other  at  the  beginning  of  stall 
(Time  E2).  Spatial  fast  Fourier  transforms  (FFTs)  of  the  pressure 
profile  are  shown  in  Fig.  12.  The  most  dominant  harmonic  is  the 
36th  mode  (its  magnitude  is  at  least  one  order  of  magnitude  larger 
than  those  of  others),  which  corresponds  to  the  number  of  rotor 
blades.  This  disturbance  is  caused  by  the  periodic  flow  features 
between  blades  and  thus  has  the  length  scale  associated  with  the 
rotor  blade  pitch.  Not  shown  here  is  the  strength  of  the  46th  mode, 
which  is  negligible.  This  indicates  that  the  influence  of  the  stator 
on  the  potential  field  at  rotor  leading  edge  can  be  ignored.  There¬ 
fore  we  will  focus  the  study  only  on  the  rotor. 

Modal  disturbances,  if  any,  should  exist  in  the  low  harmonics 
with  longer  wavelengths  that  are  on  the  order  of  the  rotor  circum¬ 
ference.  At  Phase  1  (D),  only  the  36th  mode  prevails.  This  indi¬ 
cates  that  when  away  from  stall,  the  only  disturbance  to  the  oth¬ 
erwise  uniform  annulus  flow  is  that  from  the  periodic  rotor 
passing  and  there  is  no  modal  activity.  At  Phase  2  (Time  El),  the 
early  phase  of  the  rotating  stall,  there  exists  an  identifiable  1st 
mode  along  with  the  36th  mode.  There  are  also  weaker  but  iden¬ 
tifiable  modes  between  the  1st  and  36th  modes.  As  the  simulation 
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Fig.  9  Entropy  of  Stage  35  during  stall  inception,  f=  1 07":  Left 
half  annulus  and  right  half  annulus,  16%  chord  upstream  of 
rotor  leading  edge 


proceeds,  bringing  the  compressor  even  closer  to  stall  at  Phase  3 
(Time  E2),  the  magnitude  of  the  first  mode  increases  even  more. 
The  magnitudes  of  the  modes  in  between  remain  relatively  similar 
to  those  in  the  previous  two  phases  but  fluctuate  over  the  spec¬ 
trum.  The  only  mode  that  grows  persistently  during  the  prestall 
phase  is  the  first  mode.  This  disturbance  rotates  at  100%  rotor 
speed.  This  observation  is  the  same  as  that  of  throttle  Setting  F, 
where  the  speed  of  rotating  mode  is  indicated  in  the  pressure 
history  of  Fig.  7.  This  finding  is  similar  to  the  experiment  (see  Fig. 
2)  in  that  the  one  rotor  circumference  disturbance,  which  sits  on 
top  of  the  surge  mode  (not  modeled  in  this  simulation),  starts  to 
rotate  at  80%  rotor  speed  right  before  stall.  It  suggests  perhaps 
that  the  simulation  is  adequately  modeling  the  one  rotor  circum¬ 
ference  modal  disturbance  before  the  compressor  enters  rotating 
stall. 

A  possible  explanation  of  the  origin  of  the  first  mode  can  be 
derived  from  the  following  observations:  (1)  While  the  grid  in  our 
simulation  does  not  include  the  complete  ducting  system,  it  does 
include  the  limited  ducting  that  contains  the  rotor  and  the  stator. 
(2)  The  speed  of  the  rotating  mode  is  100%  rotor  speed;  in  other 
words,  the  first  mode  is  fixed  to  the  rotor.  This  makes  it  likely  that 
the  disturbance  is  associated  with  the  rotor.  The  two  observations 
suggest  that  the  first  mode  may  originate  from  a  compound  effect 
of  the  rotor  blades  and  the  limited  volume  effect  present  in  the 
ducting  of  the  current  grid. 

Development  of  Mode  and  Spike.  To  better  understand  the 
development  of  the  first  mode  and  the  spike  disturbance,  we  focus 
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on  the  evolution  of  the  flow  during  the  transient  right  before  stall. 
Since  stall  is  likely  to  be  associated  with  the  reverse  flow  (nega¬ 
tive  axial  flow),  we  will  monitor  the  development  of  the  negative 
axial  velocity  of  the  flow  during  the  transition.  Figures 
13(a)-13(c)  illustrate  the  axial  velocity  of  the  stage  at  three  radial 
cut  planes  at  the  three  time  instants  El,  E2,  and  E3  during  the 
transient.  The  three  cut  planes  are  located  at  4  clearance  heights 
(4 s)  below  casing  (Cut  plane  1),  1.5  clearance  heights  (1.5e)  be¬ 
low  casing  (Cut  plane  2),  and  0.5  clearance  heights  (0.5e)  below 
casing  (Cut  plane  3).  Figure  13(d)  illustrates  the  front  views  of 
Cut  plane  3  at  the  three  time  instants.  The  contours  of  the  axial 
velocity  have  been  truncated  so  that  only  the  negative  velocity  has 
a  color.  All  positive  velocities  are  colored  white. 

At  Time  El,  Fig.  13(a),  about  2.5  revolutions  before  break¬ 
down,  the  reverse  flow  is  seen  in  all  rotor  passages  with  a  radial 
coverage  of  at  least  four  clearance  heights  below  the  casing,  as 
indicated  by  the  reverse  flow  being  visible  on  all  cut  planes.  The 
rotor  reverse  flow  is  formed  by  the  tip  clearance  vortex  flow, 
which  will  be  explained  in  detail  later.  The  boundary  and  magni¬ 
tude  of  the  reverse  flow  region  depicts  the  boundary  and  the 
strength  of  the  tip  vortex.  The  patterns  of  the  reverse  flow  are  not 
the  same  in  every  passage,  with  the  reverse  flow  zone  larger  on 
the  back  side  of  the  view  (210  deg).  This  can  be  readily  seen  on 
Cut  planes  1  and  2.  The  overall  shape  of  rotor  reverse  flow  cor- 


Fig.  11  Instantaneous  pressure  profile  at  rotor  tip 
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Fig.  12  FFTof  Fig.  11 


responds  well  with  the  harmonics  of  Fig.  12  in  that  the  1st  mode 
coexists  with  the  36th  mode.  The  speed  of  the  1st  mode  is  the 
same  as  that  of  the  36th  mode,  both  rotate  at  the  rotor  speed.  The 
flow  at  this  phase  has  already  shown  the  instability. 

At  Time  E2,  by  comparing  the  front  view  of  Cut  plane  3  at 
Times  El  and  E2  [Fig.  13(d)],  we  see  that  the  front  boundary  of 
the  reverse  flow  moves  forward  and  becomes  perpendicular  to  the 
axial  direction.  The  magnitude  of  the  negative  velocity  at  Time  E2 
is  also  higher  than  that  at  Time  El  as  is  shown  by  more  green 
spots  at  Time  E2.  This  is  consistent  with  the  short-length  stall 
inception  mechanism  described  by  Hoying  et  al.  [7]  that  spike 
stall  initiates  when  the  tip  clearance  vortex  trajectories  align  per¬ 
pendicularly  to  the  axial  direction.  The  larger  reverse  flow  zone 
now  rotates  to  the  front  side  of  the  view  (0  deg),  Fig.  13 (/?),  better 
seen  on  Cut  planes  1  and  2.  Also  seen  on  the  two  cut  planes  is  the 
appearance  of  reverse  flow  in  the  stator.  Different  from  the  reverse 
flow  in  the  rotor,  which  is  caused  by  tip  clearance  vortex,  the 
stator  reverse  flow  is  caused  by  suction  surface  flow  separation. 

At  Time  E3,  1.1  revolutions  after  Time  E2,  the  flow  formally 
enters  the  rotating  stall.  A  zone  of  extended  reverse  flow  covering 
five  rotor  passages  can  now  be  seen  in  the  front  view  of  Cut  plane 
3  (330  deg)  in  Figs.  13(c)  and  13(d).  This  zone  sits  in  the  middle 
of  the  less  extensive  reverse  flow  characterized  by  the  green  col¬ 
ored  tip  clearance  flow  on  Cut  plane  3.  This  less  extensive  reverse 
flow  correlates  with  the  first  mode.  The  five-passage  reverse  flow 
zone  initiates  a  spikelike  disturbance  on  top  of  the  first  mode  and 
finally  leads  to  the  rotating  stall. 

Further  analysis  of  the  simulation  also  indicates  the  correlation 
between  the  first  mode  disturbance  of  axial  momentum  and  the 
appearance  of  spike.  The  circumferential  profile  of  the  axial  mo¬ 
mentum  near  the  rotor  tip  leading  edge  is  low-pass  filtered  to 
remove  high-frequency  noise  so  only  the  first  mode  stays.  Figure 
14  shows  the  results  at  Times  El,  E2,  and  E3.  The  magnitude  of 
the  first  mode  disturbance  grows  from  El  to  E3  as  it  rotates 
around  the  annulus.  The  locations  of  the  trough  of  the  waves 
(marked  in  Fig.  14)  are  where  the  flow  has  the  lowest  axial  mo¬ 
mentum.  These  locations  match  the  regions  of  the  larger  reverse 
flows  in  all  three  cut  planes  mentioned  previously.  This  is  consis¬ 
tent  with  the  observation  of  Day  [5]  that  spikes  can  grow  out  of 
the  trough  of  the  axial  velocity  modal  waves  if  critical  incidence 
is  reached. 

As  we  pointed  out  before,  it  is  difficult  to  discern  the  actual 
transition  from  the  rotating  modal  wave  to  the  start  of  the  rotating 
stall  in  Fig.  2.  The  simulation  provides  the  details  to  complement 
the  experiment.  It  is  demonstrated  here  that  a  spike  of  five  pas¬ 
sages  long  does  occur.  It  evolves  smoothly  out  of  the  rotating 
mode  and  eventually  leads  to  the  breakdown  (rotating  stall). 
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Fig.  13  Negative  axial  velocity  at  three  cut  planes,  three  time  instants:  (a)  Time  El  (2.57  before  rotating  stall),  ( b ) 
Time  E2  (beginning  of  rotating  stall),  and  (c)  Time  E3  (during  rotating  stall);  Left  to  right,  cut  plane  1  (4e),  cut  plane 
2  (1.5e),  cut  plane  3  (0.5s:).  ( d)  Left  to  right,  front  views  of  cut-plane  3  (0.5s)  at  three  time  instants,  El,  E2,  and  E3, 
negative  axial  velocity. 


Interaction  of  Shock  and  Tip  Clearance  Flow.  Interaction 
between  the  tip  clearance  flow  and  shock  in  transonic  compressor 
has  been  cited  as  one  of  the  initiators  of  spike  disturbance.  As 


noted  by  Day  [5],  short-length  spikes  can  be  viewed  as  embryonic 
stall  cells  with  flow  breakdown  in  local  regions;  it  is  a  3D  phe¬ 
nomenon.  Computational  simulations  using  a  single  blade  passage 
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Fig.  14  First  mode  disturbances  of  axial  momentum  at  rotor 
tip 


to  study  the  effect  of  the  tip  clearance  flow  on  the  onset  of  stall 
were  reported  by  Adamczyk  et  al.  [26].  One  generally  agreed 
conclusion  is  that  the  tip  clearance  vortex  moves  forward  from 
stable  to  stall  condition.  Using  vortex  kinematics  arguments,  Hoy- 
ing  et  al.  [7]  developed  a  criterion  for  stall  inception  when  the 
trajectory  of  the  tip  clearance  vortex  becomes  perpendicular  to  the 
axial  direction.  Yo  et  al.  [8]  also  demonstrated  that  stall  inception 
is  most  likely  accompanied  by  the  forward  spillage  of  tip  clear¬ 
ance  vortex.  One  might  ask  the  following:  How  does  the  shock 
interact  with  the  tip  clearance  flow?  Does  the  shock  directly  cause 
flow  separation  that  may  contribute  to  the  stall  inception? 

To  answer  these  questions,  we  take  a  closer  look  at  the  flow 
within  the  rotor  passage.  A  comparison  of  pressure  contours  near 
the  casing  between  the  stable  and  near- stall  (Time  El)  conditions 
is  made  in  Fig.  15.  In  the  stable  condition,  the  pressure  fields  in 
the  two  passages  shown  are  almost  identical.  There  exists  a  dis¬ 
tinguishable  attached  oblique  passage  shock,  which  originates 
from  the  leading  edge  of  the  rotor  blade  and  extends  across  the 
passage  to  80%  chord  on  the  suction  surface  of  the  adjacent  blade. 
At  Time  El,  the  pressure  fields  between  passages  are  no  longer 
identical.  Even  though  the  flow  structure  is  not  identical  among 


A 


(a)  (b) 

Fig.  15  Pressure  contour  near  casing:  (a)  Stable  condition 
and  ( b )  near-stall  condition  (Time  El) 
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Fig.  16  Streamlines  forming  the  tip  clearance  vortex:  (a) 
Stable  condition  and  ( b )  near-stall  condition  (Time  El) 


passages,  there  is  still  similarity  between  them  before  the  flow 
breakdown.  We  will  discuss  the  common  features  that  are  generic 
for  the  rotor  passages  before  stall.  Examining  the  pressure  field  at 
Time  El,  we  see  that  the  single  attached  oblique  shock  now 
changes  to  a  more  complicated  dual-shock  system:  A  strong  de¬ 
tached  shock  followed  by  a  weaker  secondary  shock.  The  strong 
detached  shock,  together  with  the  stagnant  flow  around  the  lead¬ 
ing  edge,  creates  a  high  pressure  region  near  the  leading  edge  of 
the  pressure  surface.  The  pressure  outside  this  area  is  relatively 
low  until  aft  of  the  secondary  shock  where  the  pressure  level  rises 
again  to  fill  uniformly  between  the  blades. 

Streamlines  forming  the  tip  clearance  vortex  for  both  stable  and 
near-stall  (Time  El)  conditions  are  shown  in  Fig.  16.  The  pressure 
contour  near  the  casing  is  overlaid  in  the  background  to  help  the 
discussion.  First,  take  a  look  at  the  stable  condition.  Only  the 
forward  15%  portion  of  the  tip  clearance  flow  emerges  to  form  the 
core  of  the  tip  clearance  vortex.  The  trajectory  of  the  tip  clearance 
vortex  is  not  perpendicular  to  the  axial  direction.  The  oblique 
shock  intersecting  the  vortex  only  slightly  alters  the  shape  of  the 
vortex.  The  vortex  moves  away  from  the  pressure  surface  of  the 
adjacent  blade  without  impinging  it. 

At  the  near- stall  condition  (Time  El),  Fig.  16 (/?),  the  core  vor¬ 
tex,  again  formed  by  the  forward  15%  of  the  tip  clearance  flow, 
now  moves  forward.  The  center  portion  of  the  vortex  trajectory 
becomes  perpendicular  to  the  axial  direction.  This  trend  follows 
well  with  the  stall  inception  criterion  of  Hoying  et  al.  [7].  How¬ 
ever,  the  shape  of  the  vortex  trajectory  is  first  pushed  backward  by 
the  high  pressure  near  the  pressure  surface  leading  edge  and  then 
forward  by  the  high  pressure  behind  the  secondary  shock  as  it 
proceeds  downstream.  This  distortion  of  the  vortex  trajectory  is 
caused  by  the  dual-shock  system.  The  criterion  of  Hoying  et  al. 
[7],  developed  for  low-speed  machines,  still  holds  for  areas  not 
affected  by  the  dual  shock.  It  is  unclear  whether  or  not  the  dis¬ 
torted  vortex  can  accelerate  the  flow  breakdown,  but  it  is  clear  that 
the  distortion  occurs  as  the  flow  is  close  to  stall  with  the  existence 
of  the  dual  shock. 

Flow  is  likely  to  separate  due  to  the  shock-boundary  layer  in¬ 
teraction,  in  which  the  low  momentum  flow  in  the  boundary  layer 
is  subjected  to  the  sudden  pressure  rise  across  the  shock.  Thus,  it 
is  of  interest  to  examine  whether  shock-induced  flow  separations 
occur  in  the  simulation.  There  are  two  possible  locations  for  flow 
separation  in  the  rotor  passage:  One  is  on  the  suction  surface  and 
the  other  is  on  the  casing  surface.  By  examining  the  streamlines  in 
Fig.  16,  we  rule  out  the  suction  surface  separation.  To  search  for 
the  casing  separation,  the  velocity  vector  in  the  meridional  view  is 
examined.  Figure  17(a)  shows  such  a  plane  30%  pitch  off  the 
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Fig.  17  Velocity  vector  on  a  meridional  plane  30%  pitch  off 
pressure  surface:  (a)  Stable  condition  and  ( b )  near-stall  condi¬ 
tion  (Time  El) 


pressure  surface.  There  is  no  observable  reverse  flow  in  this  view. 
This  indicates  that  the  vortex  moves  straight  downstream  and  that 
the  flow  remains  attached  to  the  casing.  At  Time  El,  the  reverse 
flow  associated  with  the  distorted  vortex  can  be  seen  in  Fig.  11(b). 
The  velocity  vector  shows  that  the  flow  reversal  occurs  beneath 
the  casing  and  is  caused  by  the  organized  spiraling  fluids  of  the 
distorted  vortex.  There  is  no  indication  of  casing  separation  at  this 
instant.  The  cross  section  of  the  vortex,  however,  is  larger  than  the 
stable  condition.  This  also  shows  that  the  negative  axial  velocity 
of  the  rotor  flow  shown  in  Fig.  13  is  caused  by  the  spiraling  fluids 
of  the  tip  clearance  vortex  and  has  no  association  with  casing  or 
blade  separation.  With  the  above  observations,  we  conclude  that 
there  is  no  shock-induced  separation  in  this  during  the  stall 
inception. 

As  the  compressor  proceeds  more  into  stall  between  Times  E2 
and  E3,  the  flow  breaks  down,  i.e.,  the  mass  flow  rate  through  the 
compressor  is  declines.  The  reverse  flow  in  a  few  passages  com¬ 
bines  to  form  the  large  blockage  as  shown  in  Figs.  13 (b)  and 
13(c),  which  in  turn  expels  flow  to  adjacent  blades.  The  blockage 
creates  the  spike-type  disturbance  over  a  few  passages.  The  ex¬ 
pelled  flow  entering  the  adjacent  blade  passages  brings  in  extra 
axial  momentum  to  the  casing  areas  of  these  passages.  The  addi¬ 
tional  axial  momentum  behaves  like  tip  injection,  which  reduces 
rotor  incidence  levels,  relative  to  Time  El.  The  vortex  shape  and 
shock  structure  of  two  passages  at  Time  E2  are  compared  in  Fig. 
18.  Figure  18(a)  is  a  passage  where  the  spike  occurs.  The  shape  of 
the  vortex  trajectory  and  shock  structure  is  similar  to  those  at 
Time  El,  Fig.  16(b),  only  that  the  vortex  continues  to  grow  and 
the  high  pressure  behind  the  detached  shock  pushes  the  vortex 
backward  so  the  trajectory  moves  away  from  being  perpendicular 
to  the  axial  direction.  Spillage  of  tip  clearance  forward  of  the 
leading  edge  can  be  seen  in  Fig.  19(a),  where  it  is  colored  by  the 
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(a)  (b) 

Fig.  18  Tip  clearance  vortex  at  Time  E2:  (a)  Deteriorating  pas¬ 
sage  and  (b)  improving  passage 


pink  stream  lines.  In  the  passage  where  the  flow  condition  im¬ 
proved,  as  shown  in  Fig.  18(Z?),  the  dual-shock  system  returns  to  a 
single  shock  and  the  vortex  trajectory  is  not  distorted.  There  is  no 
spillage  of  the  tip  clearance  flow  forward  of  the  leading  edge,  Fig. 
19  (b). 

Simulation  Summary 

1.  The  simulation  results  of  the  high-speed  axial  compressor 
Stage  35  provide  evidence  that  the  current  numerical  procedure 
may  be  adequate  for  predicting  the  onset  of  flow  instabilities  and 
their  subsequent  growth  into  a  fully  developed  rotating  stall.  The 
simulation  results  clearly  show  a  full-annulus  first  mode  distur¬ 
bance  starting  to  travel  at  100%  rotor  speed  at  the  peak  of  the 
pressure  rise  transitioning  to  a  spike  disturbance  traveling  at  84% 


Fig.  19  Tip  clearance  flow  showing  spillage  over  rotor  leading 
edge  at  Time  E2  (a)  spillage  in  deteriorating  passage  and  (b)  no 
spillage  in  improving  passage 
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rotor  speed.  This  spike  disturbance  then  splits  into  multiple  stall 
cells  while  their  speeds  reduce.  Finally  the  multiple  stall  cells 
coalesce  into  a  large  single  stall  cell  that  rotates  at  43%  of  rotor 
speed  when  fully  developed.  This  simulation  models  the  natural 
growth  of  the  stall  process  without  the  use  of  an  artificial  inlet 
distortion. 

2.  As  the  simulation  is  throttled  to  stall,  small  amplitude  long- 
length  waves  of  1st  mode  and  high-order  modes  between  the  1st 
and  the  36th  can  be  discerned.  The  magnitudes  of  the  high-order 
modes  fluctuate  but  do  not  grow  much  during  the  transition  to 
stall.  The  first  mode  wave,  however,  persistently  grows  before 
breakdown.  A  short-length  disturbance  covering  five  rotor  pas¬ 
sages  bursts  into  a  spike  like  disturbance  during  the  breakdown. 
This  spike  disturbance  then  grows  into  a  rotating  stall. 

3.  The  computed  first  mode  disturbance  has  a  length  scale 
closely  resembling  the  one  rotor  circumference  disturbance  prior 
to  stall  in  the  experiment.  The  computed  speed  of  the  first  mode 
disturbance  (100%  rotor  speed)  is  also  similar  to  the  measured 
speed  (80%  rotor  speed)  of  the  first  mode  disturbance  from  the 
experiment,  estimated  from  the  results  of  Fig.  2.  The  difference  is 
within  the  resolution  of  the  picture  quality  of  the  published  data. 

4.  The  first  mode  of  the  axial  momentum  correlates  well  with 
the  appearance  of  the  growing  reverse  flow.  The  trough  of  the 
axial  momentum  matches  the  region  of  the  larger  reverse  flow  that 
leads  to  the  birth  of  the  initial  spike. 

5.  The  interaction  of  the  shock  and  tip  clearance  vortex  is  ex¬ 
amined.  As  the  compressor  approaches  stall,  the  trajectory  of  the 
tip  clearance  vortex  moves  forward  to  become  perpendicular  to 
the  axial  direction.  The  shock  system,  in  the  meantime,  also 
moves  forward,  changing  from  an  attached  oblique  single  shock  to 
a  dual  shock.  The  pressure  field  of  the  dual-shock  system  distorts 
the  trajectory  of  the  tip  clearance  vortex,  creating  a  large  reverse 
flow  region  in  the  distorted  vortex.  The  reverse  flow  in  the  rotor  is 
the  outward  tip  clearance  vortex  as  no  shock-induced  casing  or 
blade  separation  is  found. 

6.  During  the  flow  breakdown,  spillage  of  the  tip  clearance  flow 
forward  of  the  rotor  leading  edge  is  found  in  several  passages, 
where  the  flow  condition  is  deteriorating.  A  large  blockage  region 
is  formed  in  the  deteriorating  passages;  this  in  turn  impels  flow  to 
enter  into  adjacent  passages.  The  influx  of  fluid  momentum  in  the 
adjacent  passages  improves  the  flow  condition  there. 


Conclusions 

The  turbo  code  has  been  applied  to  the  3D  full-annulus  simu¬ 
lation  of  a  transonic  compressor  stage  NASA  Stage  35  to  investi¬ 
gate  the  flow  evolution  prior  to  stall.  This  simulation  does  not 
include  the  system  volume,  and  therefore  cannot  model  the  surge 
mode  oscillation  that  arises  from  the  test  facility  ducting  system. 
With  the  surge  mode  removed,  the  stalling  flow  first  takes  the 
form  of  a  rotating  long-length  (modal)  disturbance;  it  then 
changes  into  a  spike.  The  spike  seems  to  grow  out  of  the  long- 
length  disturbance  and  eventually  leads  to  stall.  The  spike 
emerges  when  the  tip  clearance  vortex  becomes  perpendicular  to 
the  axial  direction  and,  in  some  passages,  with  forward  spillage. 
The  passage  shock  changes  from  single  shock  to  dual  shock  when 
approaching  stall.  The  dual  shock  distorts  the  trajectory  of  the 
vortex  but  does  not  cause  flow  separation  that  directly  contributes 
to  the  stall.  The  simulation  correlates  well  with  the  experiment 
data  excluding  the  surge  mode.  The  computational  resources  in¬ 
volved  in  conducting  such  a  simulation  are  prohibitively  large  at 
the  time  and  are  not  suitable  for  routine  executions.  While  this  is 
certainly  an  issue,  the  key  point  is  that  the  simulations  of  this  type 

Journal  of  Turbomachinery 


generate  qualitatively  and  quantitatively  accurate  flow  field  data 
that  can  support  our  understanding  and  inferences  of  the  unsteady 
compressor  flow  field,  including  support  and  interpretation  of  ex¬ 
perimental  work. 
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